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Abstract 

We study the asymptotic behaviour of a spinning top whose shape is spheri- 
cal, while its mass distribution has axial symmetry only, and which is subject to 
sliding friction on the plane of support (so-called tippe top). By a suitable choice 
of variables the equations of motion make explicit the conservation of Jelett's in- 
tegral (rediscovered by Leutwyler) and allow to construct explicitly all solutions 
of constant energy. The latter are the possible asymptotic states of the solutions 
with arbitrary initial conditions. Their stability or instability in the sense of Lia- 
punov is determined for all possible choices of the moments of inertia. We conclude 
with some numerical examples which illustrate our general analysis of Liapunov 
stability. 



1. Introduction 

The so-called tippe top may be modeled by a sphere whose mass distribution is axially 
symmetric, but not spherically symmetric, so that its center-of-mass does not coincide 
with its geometrical center. As described in the historical introduction to the article 
by R.J. Cohen on the subject this top's astonishing behaviour has puzzled several 
generations of physicists: Provided certain conditions on the moments of inertia I\ = I2 
and J3 are fulfilled, the rapidly spinning top will quickly tip over to a completely inverted 
position where the center-of-mass sits vertically above the geometric center. Thus, in the 
initial phase of the motion the center-of-mass is raised, with the rotational and the total 
energy decreasing to a constant value. In the same process the direction of rotation, 
with respect to a body-fixed frame, is reversed. 

A detailed analysis of the tippe top and a numerical study of its equations of motion 
was presented by Cohen |]|. In particular, this analysis established definitely an earlier 
contention by Pliskin, Braams, and Hugenholtz 0, |3|, Q]: It is the sliding frictional 
force acting at the point of contact between the top and the plane of support which is 
responsible for the inversion. 

Recently, Leutwyler showed that if that frictional force is the dominant one and if 
it is proportional to the sliding velocity then, while the individual projections L3 and 
L3 of the angular momentum onto the spatial vertical and onto the body's symmetry 
axis, respectively, decrease by dissipation, the specific linear combination A = L 3 — aL 3 , 
(with a denoting the distance of the center-of-mass from the sphere's center, in units of 
its radius), remains constant in time [||]. This conservation law is of a purely geometric 
nature and does not depend on the specific form of the frictional force, as a function 
of velocity, (for a proof using the geometry of the system only, see 0). Using this 
constraint, Leutwyler showed, by a simple energy consideration, that the inverted state 
of the spinning top, indeed, has rotational energy lower than the non-inverted state. 

In fact, that conservation law, called Jelett's integral was known to the experts 
in the field much earlier [H]. In these references as well as in the dynamics of gyro- 
scopes was studied. Synge, in particular, studied the stability of an asymmetric tippe 
top spinning in inverted position ||. 

To the best of our knowledge none of these articles addressed the long-term orbital 
stability of this system, in the sense of a Liapunov analysis. The questions to be asked 
being: Which of the possible inverted positions (whose nature depends on I\ and I3) 
is asymptotically stable? Which initial configurations are driven into these asymptotic 
states, and how do they move towards them? In this work we present a rigorous and 
complete analysis of the (symmetric) tippe top. We study its long-term behaviour by 
means of a Liapunov function, making use of the conservation law, and thereby give 
a complete description of the asymptotic solutions and of their stability or instability. 
The conservation law, clearly, allows to reduce the number of variables, as compared 
to Cohen's analysis. Furthermore, by a suitable choice of these variables the equations 
of motion are simplified further. Repeating then the numerical study, the pattern of 



2 



explicit solutions is rendered considerably more transparent. 

In sec. 2 we define the system that we study, in more precise terms, list our assump- 
tions, and state the conservation law. In sect. 3 we write down the equations of motion 
and derive all solutions with constant energy, i.e. the asymptotic solutions towards which 
the spinning top will tend if they are found to be stable. Sect. 4 addresses this stability 
analysis for any initial condition. This section contains our main results. In sect. 5 we 
give some examples of our own numerical study and summarize our results. 

2. Definitions, assumptions, and conservation law 

The top is taken to be a sphere with unit radius, r = 1, and an axially symmetric 
mass distribution. Its symmetry axis is taken to be the body-fixed 3-axis so that for 
the moments of inertia I± — I2 ^ I3. The center-of-mass S has the distance a from the 
sphere's center M, with < a < 1. A is the point of contact with the plane of support 
as sketched in fig. 1. 

Generally speaking, there are three different types of motion which are possible: (i) 
the top rotates about a vertical axis through a fixed point on the plane. In this state of 
motion that we call rotational below, only rotational friction is active; (ii) the spinning 
top rolls over the plane without sliding. For this motion that we shall call tumbling below, 
only rolling friction (and, possibly, rotational friction) is active; (iii) more complicated 
spinning whereby the top slides over the plane of support and, hence, is subject to sliding 
friction. Among the type (iii) solutions we consider only those for which the top is in 
permanent contact with the plane, that is, we do not consider hopping states of motion. 

To the extent that rolling friction and rotational friction can be neglected as compared 
to sliding friction, solutions of type (i) or (ii) are asymptotic solutions, with constant 
energy, which may or may not be stable. We shall assume, indeed, that rolling friction 
as well as rotational friction are absent. The problem to be solved then is twofold: to 
classify all solutions with constant energy and, by means of a Liapunov analysis, to 
study the long-term behaviour of dissipative solutions with general initial conditions. 

Sliding friction slows down the projection L 3 of the angular momentum onto the 
vertical direction, i.e. the laboratory 3-axis, as well as its projection L 3 onto the top's, 
body-fixed, 3-axis, through torques R and R, respectively, due to the frictional force 
that acts at the point of contact A, viz. 



The (sliding) velocity components of the instantaneous point of support A which are 
due to infinitesimal rotations about the 3-axis and about the 3-axis, respectively, have 
the same direction in the plane of support and are perpendicular to the plane spanned 
by the 3- and 3-axes. The velocity component of A which is due to a change d9 of the 
angle between the 3- and the 3-axis, on the other hand, lies in that plane and, hence, 
is perpendicular to the former two components. Obviously, analogous statements hold 




-R, 




R. 
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true for the components of the frictional force acting on A, and are independent of its 
explicit functional dependence. As a consequence, the torques R and R depend on the 
same component of that force and differ only by the moment arms whose lengths are 
seen to be asm9 and sin^, respectively, from fig. 1. Thus, R = aR and, from eq. (1), 
the linear combination 

A := L 3 - aL 3 (2) 

is a constant of the motion. 

If the rotational kinetic energy is large as compared to the potential energy in the 
gravitational field, the energy of the spinning top can be rewritten in terms of A as 
follows 

E « T rot = A 2 - z 2 ) + h{z - a) 2 }' 1 (3) 

where z := cos 9. This formula can be used for deriving criteria for partial or complete 
inversion |J. For instance, if the moments of inertia obey the inequalities 

(l-a)/ 3 <Ji<(l + a)/ 3 (4) 

the expression (3) assumes its smallest value in the completely inverted position z = —1. 

In our analysis below we will recover the conservation law (2) from the equations 
of motion. The criterion (4), as well as analogous criteria for other types of inverted 
motion, will appear in the stability analysis, though modified in the presence of gravity. 

3. Equations of motion and solutions with constant energy 

3.1 The equations of motion 

As customary in the theory of rigid bodies the motion of the top is described most 
conveniently using three systems of reference ||: An inertial, laboratory, system K 
whose 3-axis is taken to be the vertical; a (non-inertial) system K which is attached to 
the center-of-gravity and whose axes are parallel, at all times, to the axes of Ko; and 
a body fixed, principal-axes-system K whose 3-axis is the symmetry axis of the top. If 
eg denotes the unit vector along the top's symmetry axis with respect to K, and R(t) 
the SO (3) rotation matrix which connects K and K, that same unit vector is expressed 
with respect to K as follows 

fj = R(t)e- 3 . (5) 
Again with respect to K the inertia tensor is given by 

i(t)=h{t+ I -^\m)(m\}, (6) 

in an obvious notation for the dyadic constructed from the vector r), (| )), and its 
transposed (( |). The inverse of the inertia tensor is seen to be 
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The angular velocity u>(t) is defined by the formula 



ad u>(t) = w(f)x = R(t)R T (t) = -R(i)R (t) , (8) 

while the angular momentum is L(t) = l(t)u)(t). By eq. (7) u> is expressed in terms of 
L as follows: 

u(t) = j i {L(t)-I^(r,\L)r>}. (9) 

The variables r) and L are shown in fig. 2, for the same position of the top as in fig. 1. 
Let s(t) be the coordinate vector of S, v the velocity of the momentaneous point of 
support A, F the external force in the laboratory system K , and let N denote the 
external torque in the sytem K. By the axial symmetry of the top F and N depend 
only on (r), L, s). The equations of motion read 

d „ . 1 „ 

— rj = u x 77 = — L x rj 

= N(fi,L,s) 
ms = F(fi,L,i), (10) 

with m denoting the total mass of the top. In the first of these equations we have made 
use of eq. (9). As we require the top to be in contact with the plane of support at all 
times, the component S3 of s is not an independent coordinate. Clearly, the requirement 
is that the 3-component of the coordinate vector of the momentaneous point of contact 
A must be zero at all times. With a denoting the the vector that joins the center-of-mass 
S to the point of contact A (cf. fig. 2) 

a = aff — e 3 (11) 

this implies 

(e 3 I s + a) = s 3 + a(e | 77) - 1 = 0. (12) 
With v and defined by 

v := s + x a , := u> x a , 

and making use of eqs. (9) and (11), there follows from eq. (12) 

v 3 = (e 3 I s + u x a) = s 3 + ^ (e 3 | L x 77) = . 

This result shows that s 3 , the 3-component of the center-of-mass' velocity, is not an 
independent variable but is a function of L and t), i.e. s 3 = s 3 (t), L). Therefore, the 
third equation of the system (10) must be replaced by msi^ = Proj^-F 1 , where Proj 12 
denotes the projection onto the plane of support. 
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It remains to derive explicit expressions for the external force F and the external 
torque N . The external force is the sum of the gravitational force F g = —mge 3 , the 
normal force F n describing the action of the plane of support at A, F n = g n e 3 , and 
the frictional force Ff whose direction is opposite to the velocity of A in the plane, 
Ff — —gfV, with g n and gj positive semi-definite functions, F = F g + F n + Ff. The 
force exerted in the point A, in turn, is the sum of F n and Ff, Fa = F n + Ff, so that 
the torque N is given by 

N = axF A = (afi-e 3 )x(g n e 3 -g f v). (13) 

The final form of the equations of motion is then 

d „ 1 
—77 = —L x 77 

at i\ 

= (afj - e 3 ) x (g n e 3 - g f v) 
ms 1>2 = -gfb. (14) 

Before we move on we note that the conservation of the quantity A, eq. (2), follows from 
the equations of motion (14). We have X(ff, L) = L 3 — aL 3 = —(a \ L) and 

-\=(a\L) + (a\N) = 0. (15) 

The first term is independent of the force of friction and, hence, vanishes because for a 
force-free, axially symmetric top both L 3 and L 3 are conserved. Alternatively, this may 
also be seen from da/dt = adf)/dt and the first of eqs. (14). The second term vanishes 
because the torque N = a x F A is perpendicular to a. As stated in the introduction this 
result is independent of the explicit functional dependence of the frictional force. The 
way we have chosen the independent variables the conservation law is already encoded in 
the equations of motion (14)Q Therefore, the seven variables that appear in the system 
of equations (14) form an optimal set of independent variables. It will become clear, 
furthermore, that this set is optimal in the sense of being well adapted to the problem 
we are studying. 

In our numerical analysis below we shall assume the frictional force to be proportional 
to the normal force, i.e. 

9f = V9n , 

with /i a constant, positive coefficient of friction jj]]. The coefficient g n , that is the 
magnitude of the normal force, is calculated from the orbital derivative of S3 and from 
Newton's law. It is found to be a function of the seven independent variables, g n = 
g n (f), L, si )2 ). The result is given in equation (A.l) of the appendix ]TD . 



1 Indeed, Cohen's analysis Q leads to 10 first-order differential equations while the system (14) 
involves only 9 equations of first order in time. 
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Finally, we note that the frictional force, if it is taken to be proportional to v/ \\v\\, 
is undefined for v = 0. As the asymptotic states of the top involve configurations where 
v vanishes, we replace the expression for the frictional force by a functional form that is 
continuous in v = and vanishes at that point. This is achieved most easily by replacing 
the unit vector v/ \\v\\ by 

v = h(\\vr v 



,i v ii 

where the function h(x) is chosen such that it fulfills the conditions 

h > , h(x) = •<=>- x — 0, | h(x) — 1 |< 5 for all x > e for given e, 5 . 

An example for such a function is h(x) = tanh(iVx) with iV a sufficiently large positive 
integer. 

Any solution {ff, L, Si^} of the system (14) for which the coefficient g n (fj, L, Si^) is 
positive, is physically admissible. Therefore, the domain of definition for the equations 
of motion is 

n = g-\]0, oo [) C S 2 x 1R 3 x 1R 2 . 

On this domain the equations of motion are real analytic. 

3.2 Solutions of constant energy 

Recall that we assume sliding friction to be the only frictional force present, rotational 
and rolling friction being neglected in our analysis of the tippe top. Clearly, any solution 
of the system (14) for which the sliding velocity v of the point of contact vanishes 
at all times, must have constant energy. The converse is also true: any solution with 
constant total energy has the property v = 0. Indeed, as we shall confirm in sect. 4, 
E = —fj,g n h(\\ v ||) || v || . Therefore, all solutions with constant energy satisfy a system 
of differential equations which follows from (14) by setting v — 0, viz. 

d „ 1 

dT 7 = h Lxr > 

d _ 

— L = ag n r] x e 3 

77l«i )2 = 0, (16) 

where the coefficient g n is given by eq. (A.l) with v = 0, 

* = *<fi,L)=rv l+ ^\-™%® . (17) 

1 + ma z (l — Tjz) /Ii 

Note that the velocity v refers to the laboratory system K , i.e. v = v 1>2 = si, 2+^1*2 — 0> 
and, from the third equation of the system (16), v^ 2 — const., referring to the 
system K which is attached to the center-of-mass. Therefore, the third equation (16) 
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(s) 

is equivalent to the condition v\ 2 — const, or, for its orbital derivative with respect to 

eqs. (16), Vil = 0, so that we may as well study solutions of the first two equations to 
which we add that subsidiary condition 

v$(fi,L) = 0. (18) 

Before giving the explicit form of the solutions with constant energy we collect their 
properties in the following 

Proposition 1: // rotational and rolling friction are absent all spinning solutions of the 
tippe top with constant total energy are characterized by the properties 

(i) the projections of the angular momentum L onto the vertical and onto the top's sym- 
metry axis are conserved, L 3 = (e 3 | L) = const., L 3 = (fj \ L) = const., 

(ii) the square of the angular momentum is conserved, L 2 = const., and so is the pro- 
jection of fj onto the vertical, r] 3 = (£3 | 77) = const., 

(Hi) at all times e 3 , fj, and L lie in a plane, 

(iv) the center- of-mass stays fixed in space, s = 0. 

Proof: (i) is a well-known result for the children's top in a gravitational field || and 
follows from the first two equations of motion (16). To prove (ii) and (iii) we must cal- 
culate v( s \ the velocity of the point A with respect to the system of reference K, and 
make use of the condition (18). 

With to as given by eq. (9) and a as given by eq. (11), one has 

fW ..... 1 r , .. . , Wi- 



v (s> = uj x a = — < aL x fj — L x e 3 H L 3 r) x e 3 V . (19) 

ill I3 ) 

From this expression one calculates the orbital derivative of v^ s \ by means of the equa- 
tions of motion (16). As the vectors e 3 x i) and e 3 x (e 3 x f)) = 77363 — r) both are in the 
(l,2)-plane, the condition (18) requires their scalar products with i)^ to vanish. This 
leads to the conditions 

Lsfa-^^T^P = 0, 
a(l-r] 2 3 ) (L 2 - g n h{\ - arj 3 )) + 
L 3 l [ a( V3 L 3 -L 3 )- I -^^(L 3 - V3 L 3 )} = 0, (20) 

where P stands for the scalar product P := (e 3 x L \ fj) (and its cyclic permutations). 
By the equations of motion the orbital derivatives of r/ 3 and of L 2 are also proportional 
to P, dr] 3 /dt = P/h, dL 2 /dt = 2ag n P. Therefore, upon multiplication of the second 
equation by P, the above conditions can be rewritten as follows 

a j-Vs) Tt Vs =0, 



8 



Igja^La-Xs)-^^^-^)}^! = 0. (21) 

As long as L3 7^ the first eq. (21) implies 773 = const., hence 773 = and P = 0. Then 
we have also d/dt(L 2 ) = 0, thus proving (ii) and (hi). The case L 3 = is a little more 
complicated: If L3 = and if r] 2 7^ 1, the second equation (21) implies that the product 

L 2 (l - ar] 3 ) 2 = const. (22) 

is a constant. The second equation (20), in turn, reduces to L 2 = g n Ii(l — ar) 3 ). Inserting 
here the expression (17) for g n gives the relation 

, 9/ ma ma 2 . 
mgh{\ - arfy) - L 2 (l - — m + — — ) = . 

Finally, replacing L 2 by means of eq. (22) we obtain a cubic equation for 773 with constant 
coefficients. This equation has at least one real solution which is a constant. 
This proves that 773 is constant in all casesQ and, thus, that the quantity P vanishes, i.e. 
that £3, L, and r) lie indeed in a plane. The velocity v^ s \ eq. (19), lies in the (l,2)-plane 
and, by eq. (12), S3 = 0. 

The last part (iv) follows from the explicit solutions that we give next. These solutions 
pertain to the following classes: 

(A): 773 = 1, (B): 773 = —1, (rotating solutions). As t) = ±£3 the first eq. (16) implies 
L = A e 3 , with A = ±V?, L 3 = A , L 3 = ±A , and, according to eq. (2), A = 
A (l =Fa). Eq. (17) reduces to g n = mg = const.. One verifies that the second condition 
(20) is fulfilled and that eq. (19) yields <u (s) = 0. 

(C): — 1 < T] 3 < +1 (tumbling solutions). As the projection 773 of the unit vector r)(t) on 
the vertical is constant, its time dependence must be of the form 

f){t) = R 3 (0(t))T) (o) , with R 3 (0) = exp{0e 3 x} . 

Furthermore, according to (iii), the angular momentum has the decomposition L = 
Ae 3 + Ar)(t), where A and A are constants. Then L 3 = A + A773, L 3 = A + Arj 3 , and, 
from eq. (19), is again seen to vanish, = 0. 

The first equation of motion (16) gives = A/Ii = const., the second equation of 
motion, together with eq. (17) yields the relations 

AA = —amgl\ , g n = mg . (23) 

2 Thus, an oscillatory or rolling motion whereby the top swings or rolls about an axis perpendicular 
to the symmetry axis and parallel to the plane of support, is subject to sliding friction and does not 
have constant energy. 
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If these are inserted into the second equation (20) one obtains 



A 2 = m9 T f (24) 

For a given value A of the conserved quantity (2) the parameter r/ 3 can be expressed in 
terms of A in two different ways, viz. 

= ajmgll - I 3 A 2 ) = A 2 - AA + o?mgI x 
V3 ~ tfih-h) ~ a(A 2 + mgh) ' [ } 

The first of these follows from eq. (24), the second form follows from the expression 
A = L 3 — aL 3 = A(l — ai] 3 ) + A(i] 3 — a) and from eq. (23). From eq. (25), finally, one 
obtains the fourth-order equation for the quantity A 

J l 2 A a 4 ^3-^1,3 , / T n 2 J 1 



— a I A — A--— —A + (amgh) j- = . (26) 

Thus, in all cases the velocity vanishes and, from the condition v = s + = 
also the 1- and 2-components of the center-of-mass' velocity vanish. This proves (iv).D 

3.3 More about the tumbling motions 

As a preparation for the stability analysis we need to know how many tumbling 
solutions there are, given the two moments of inertia, for a given value of the constant 
of the motion A, eq. (2). The answer to this question is provided by the following 
statements and results. 

Tumbling solutions: (Tl) Equation (26) which is of degree < 4 in the unknown A, has 
at most two real solutions. Indeed, the function y(x) = a 4 x 4 + a 3 x 3 + ao whose derivative 
has a double zero at x — 0, has at most one extremum. As two zeroes are separated by 
at least one extremum, y can have no more than two real zeroes. 

A solution of eq. (26) is admissible only if A is real and if rj 3 , as calculated from the 
equations (25), lies in the interval (—1,1). The following special cases are particularly 
easy: 

(T2) For I\ = I 3 the solutions of eq. (26) are A± = ±\/mgI\. For i] 3 to be in the 
admissible interval, the constant A must obey the inequalities 

(l-af< ^<(l + «) 2 - (27) 

As these can hold at most for one of the two solutions, there is at most one tumbling 
solution. The derivative y'(x = A±) = —4a 2 mgli being different from zero there is then 
a neighbourhood of the value A^ = A±, and of J 3 = I\ for which there is precisely one 
real zero of eq. (26) which guarantees the condition —1 < r] 3 < +1. 
(T3) Let Ii = J 3 (l — a 2 ). In this special case eq. (26) has exactly one real solution which 
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is Ao = ((mg/i) 2 (l — a 2 )/A) 1 ^ 3 . There is at most one type of tumbling motion. 
(T4) In the limit |A| ^> y/mgli, i.e. in the limit where the rotational kinetic energy is 
large as compared to the gravitational energy, counting the tumbling solutions becomes 
particularly simple: 

(i) If ii > 73(1 + 0:), or if Ji < 73(1 — 0;), there exists a positive number cq such that for 
|A| > Co eq. (26) has one and only one solution which yields 773 with — 1 < 773 < +1. In 
other words, there is exactly one tumbling solution. 

(ii) For J 3 (l — a) < h < J 3 (l + a) there exists a positive Ci such that for all |A| > C\ 
there is no real solution of eq. (26) satisfying the subsidiary condition —1 < r] 3 < +1. 
There are no tumbling solutions. 

The proof goes as follows: We consider first the special cases (T2) and (T3) both of which 
belong to case (ii). If I\ = I3 the inequalities (27) are violated for all |A| > A^/mgli. If 
l\ = J 3 (l — a 2 ), A is as given above, and i] 3 (as calculated from the first equation (25)) 
is in the right interval only if A lies in the interval 

1 — a A l + o 

< ; < 



\/l + a y/mgh y/l — a 



Clearly, this is not the case whenever |A| ^> \fmgT\. 

We then consider the general situation where 1\ neither is equal to I3 nor to J 3 (l — a 2 ) 
Eq. (26) which we rewrite as follows 



1 ~ l^h ) A + (m9/lQ) TT^A A 5 = A ' 



for sufficiently large |A|, has two real solutions whose asymptotics is 

A / / \ 1/3 

Al ~ 1 h 2 ' As ~ ( m 9ha) 2/3 — J — - A- : 

The first of these, A 1 tends to infinity as A — > 00 and, from eq. (25), 



ife(Ai) — j^-a (28) 
_ J-i 

which is indeed in the right interval if the inequalities (i) hold true. The second solution 
tends to zero, for large values of A, while 773 (A 2 ) tends to infinity and, hence is not in 
the right interval. This proves (ii). 

4. Long term behaviour and stability 

In this section we show that the total energy of the top is a suitable Liapunov function 
and, in case of asymptotic stability, that the solutions tend towards some solution with 
constant energy. This analysis makes essential use of the conservation law (2) and of the 
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solutions with constant energy that we studied in sees. 3.2 and 3.3 above. 

4.1 The energy is a Liapunov function 

With our choice of variables the total energy of the spinning top is given by 

£(r), L, * 1>2 ) = | [*? >2 + sl(f>, L)] + ^- (L 2 - ^=^3) + mgs 3 (f,) . (29) 

The first term is the kinetic energy of the center-of-mass motion, the second is the 
rotational energy T rot = u> • L/2, with u> as given in eq, (9), the third term is the 
potential energy. The orbital derivative of E(ff, L, s 12 ) is calculated from the equations 
of motion (10) or (14), and from eq. (12) for s 3 . The result is 

^E(fj, L, s 1>2 ) = v(j), L, a 1)2 ) • F f (f}, L, s 1>2 ) = -fjig n h(\\ v \\) \\ v \\ . (30) 

We recall that v is the velocity of A with respect to K (cf. sec. 3.1), and that g n and 
v = 51,2+^1*2 depend on r), L, and s 1)2 , cf. eqs. (A.l) and (19). The coefficient of friction 
H being positive, the orbital derivative of E(t) is negative semi-definite. It vanishes if 
and only if v(f], L, Si j2 ) vanishes. The function E(t) decreases monotonically and, hence, 
is a Liapunov function. Furthermore, the equations of motion and the function E are 
real analytic. Thus E(t) is either strictly monotonous or is a constant. As expected, 
E becomes constant whenever the sliding velocity v vanishes. As E(t) is a Liapunov 
function the asymptotic states, for t — > 00, are solutions of constant energy. 

4.2 Extrema of the Liapunov function 

In determining the extrema of the Liapunov function E(t) = E(ff(t), L(t), s lt2 (t)) on 
the hypersurfaces which are defined by the condition 

\(fi,L) = A (0) = const. (31) 

we make use of the following idea: In case of asymptotic stability the system will tend 
towards one of the solutions of constant energy. In these asymptotic states vanishes 
and, from eq. (19), the angular velocity uj is proportional to the vector a. eq. (11). 
Therefore, L is proportional to la, with I as given in eq. (6). Furthermore, L is in the 
same plane as r) aned £3. For an arbitrary spinning state of the top (in which sliding 
friction is still active) we decompose the angular momentum L into a component 

L\\ = I(f))a = yL\\ 1(f)) (of, - c 3 ) (32) 

parallel to what its direction would be if E were already constant, and a component L± 
perpendicular to a, Lj_ ■ a = 0. With a extracted from eq. (32) this means, in fact, that 
L\\ and L± are orthogonal with respect to the scalar product defined by I _1 (f)), viz. 

T _ u L„ 



(L ± \r'(f ] )\L ll } = -^L ± .a = 0. 
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When the spinning state tends towards a solution of constant energy, L± will tend to 
zero, L will tend to its asymptotics L\\. 

With L — L\\ + L±, Ly being defined by eq. (32), we have 

(u\L) = | L) = (r 1 ^ | L\\) + {r 1 L ± I L ± ) 

and (e 3 | L x fj) = (e 3 | L± x 77). The total energy can be written as the sum of two 
terms 

E = E^ 1 \ m ,L ll ) + E^(f,,L ± , Sl , 2 ), (33) 
the second of which contains all terms that will vanish asymptotically 

= \{I- l L± I L ± ) + I [(s li2 I s li2 > + a 2 (f, x e 3 | J L ± ) 2 // 1 2 ] , (34) 

while the first depends on 773 and L\\ only, and is given by 

L 2 

E W = -J-G( m ) + mg(l - arj 3 ) (35) 



with 



G (V3) = 7-(a I 1 I a) = -jr(afi-e 3 \ I\afj-e 3 ) 
H I i 



= l-r& + y(m-<*) 2 . (36) 



The constant of the motion (2) is given by 



A = -(a I L) = -^(a \I\a) = -L\\G(i] 3 ) 
so that L\\ can be calculated from i] 3 and A, viz. 

LlM = -ak)' (37) 

which, by proposition 1 (ii) becomes a constant for t — > 00. 

Clearly, both E(ff, L, 8^2) and X(fj, L) are invariant under rotations R 3 (0) about the 
vertical. Therefore, E will be extremal under the subsidiary condition A = \^ = const, 
at most on the sets 

r(rj, L, si )2 ) := {(R 3 (0)^7, R 3 (0)£, R 3 (0)si, 2 ) | G [0, 2vr]} . 

The following proposition shows that determining the extrema of the function E, eq. (33), 
on the hypersurface defined by the condition A = A^, in fact, is equivalent to finding 
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the extrema of E^ l > as a function of 773 and, thus, is a one-dimensional problem^. We set 
773 = cos# and write E^(9) for E^(rj 3l L\\{r) Z: A (0) )). Then, 
Proposition 2: The following assertions are equivalent: 

(i) the function E{f], L, Si^) assumes an extremum under the condition X(f],L) = X^ 
for(f 1 ,L,s h2 )er(rt \L( \s?) l ); 

(ii) and sf'l have the values 

L<°> = L® = AW) I(f)(°))(^W - e 3 ) , = (38) 

and i/ie function E^(6) has an extremum for cos 6 = 773° . In particular, the minima of 
E^{9) correspond to minima of E on the hypersurface X(i),L) = \(°>, while maxima or 
saddle points of E^ yield saddle points of E with A = X^ . 

The proof is easy and we do not write it down here fH|. It makes use of eq. (34) which 
shows that E^ 2 > is a strictly monotonous function of || L± || and of || Si^ ||, and vanishes 
when these vectors vanish. 

Finding the extrema of the function E^(9) as defined by eq. (35), with A = A 1 - ** fixed, 
with G(r/3 = cos#) as given by eq. (36) and with L\\ = —X^/Gfa), eq. (37), is tedious 
but straightforward. We merely give the result in 
Proposition 3: Let 6 = arccosr^ '' be such that dE^(6)/d6 = 0. 

(i) If 7/3^ = ±1, E^(6) assumes a minimum (maximum) iff 

I 3 (l T a)-J 1± I^(l T a) 4 (39) 

is positive (negative). 

(ii) If —1 < 773° < +1, E^(6) assumes a minimum (maximum) iff 

Js(l - a 2 )/ [l + 3((I 3 - h)rtf> - ahf/I 2 ] - h (40) 
is negative (positive). 

On the basis of proposition 3, and taking into account the condition (38) as well as 
eq. (37) for Lu, it is plausible that the set T C on which E is extremal, with A = X^ 
fixed, coincides with the equivalence classes of solutions of constant energ 

yft That this 

is indeed so is the content of proposition 4 whose proof we again omit [K| . 
Proposition 4: The sets T(f)^°\ L^, s^ 2 ) for which and sf^ are given by eq. (38) and 
for which E^' = at 6 = arccos 773 , are precisely the trajectories with E(t) = const, 
and X(t) = X^ . In particular, the cases r]^ = +1 and —1 correspond to the classes (A) 
and (B) above (rotating solutions), respectively. The case —1 < r/f^ < +1 corresponds 



3 This generalizes Leutwyler's result who neglected gravity H. 

4 Two solutions are called equivalent if they have the same trajectory. 
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to the tumbling solutions (C). Since for a given value \^ of A there exist at most two 
tumbling solutions, E^(6) can have at most two extrema, that is, can have at most one 
minimum on the open interval < 6 < tt. 

In the first two cases, 773 = ±1, we have from sec. 3.2 

T(^°\L^,s^ 2 ) = {(±e 3 ,A e 3 ,0)} . 

In the second case = 0> an d an d are obtained from eqs. (23) - (26). 
4.3 Asymptotics of solutions and Liapunov stability 

Let $t denote the flux of the equations of motion (14), or, for E = const., (16), 
(fj(°\ L,(°\ sf'l) the initial conditions. We have shown above that if the energy is con- 
stant, E =const., and if A = A'- -' is given, there are at most four different trajectories. 
We denote the time-positive trajectories of constant energy, if they exist, by 

7± = {$ t (±e 3 ,Aoe 3 ,0) | t > 0} ; 7l = {$ 4 (r(r)«, L®, 0)) | t > o} . (41) 

The first of these are the rotating solutions, the second are the tumbling solutions with 
— 1 < 773 < +1 that we studied in sec. 3.3. 

Consider now the solution of the full equations of motion (14) pertaining to the 
arbitrary initial condition (ff^°\ L^°\ sfl) and let u denote its limit setf]. The asymptotic 
behaviour of the solution is fixed by the following proposition. 

Proposition 5: Let (ff(t), L(t), s li2 (t)) be the solution defined on the interval l ma x C IRt, 
with initial condition (f)(0) = rf \ L(0) = -L*' -', Si^O) = s^l). Assume further that 
there is a positive constant g^ such that 

g n (f](t), L(t), s 1>2 (t)) > g^ > , for all t G \ max . 

Then l max D [0, oof and there exists exactly one trajectory 7 G {7+, 7_, 71, 72} such that, 
as t — > 00, the solution tends to 7. The 00 limit set of (r)^, L^°\ s^l) is 7. 
We sketch the proof: For all t > the energy is bounded E(t) < E(t = 0). There- 
fore, 7 ( - + ' 1 = |$t(T) < -°- ) , L^°\ Big) I t > o| is bounded. By assumption the solution is con- 
tained in the set A4 = |(r),X, s ii2 ) G Q \ g n > g^} which is a closed subset of fi. 
Therefore l max C [0, 00 [ and 7'+' is relatively compact. Following standard theory of 
ordinary differential equations ]TT[ we conclude that (i) 00(7)^, L^°\ s^l) is compact, 
connected and positive invariant by $ t ; (ii) there exists a real value E^ such that 
cj(r) (0) , £ (0) , sjjj) = E-\E^); (iii) for t -> 00 the solution tends to u(f] {0 \ L {0) ,s { ° ] 2 ). 
The statements (i) and (ii) imply that this limit set is the connected closure of tra- 
jectories with E = E^ =const. and A = A^ = X(fj ( - \ L^). Knowing that any two 

5 Recall that the u limit set of a;'- -' for the flux <I>t is the set of all accumulation points of <&t(x^) as 
*-> 00, cf. e.g. 0. 
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of these trajectories have a finite distance we see that there is exactly one trajectory 
7 G { 7+ , 7-, 7i, 72} such that L®, a$) = 7 ^ 0- 

We now turn to the central topic of our investigation: the question of orbital stability 
of the spinning motions of the tippe top. The answer is contained in the following 
Theorem : 

(i) If the quantity (39) with the upper sign is positive then 7+, the non-inverted, rotating 
motion, is Liapunov stable. If it is negative, 7+ is unstable. 

(ii) If the quantity (39) with the lower sign is positive then the completely inverted, 
rotating motion, is Liapunov stable. If it is negative 7_ is unstable. 

(Hi) Let the quantity (40) be negative (positive). Then in as much as the tumbling mo- 
tion corresponding to 773 exists, ji is Liapunov stable (unstable). 

To prove this theorem let 7 be 7 = 7+ and 7 = 7_, in the cases (i) and (ii), respectively, 
and let 7 = 7^ in the case (iii). According to proposition 3 the conditions given in the 
theorem are sufficient for E^ l > to assume a minimum (maximum) for cos# = 1 and — 1, 
for (i) and (ii), respectively, or cos 6* = 773 for (iii). By proposition 2 this means that 
the total energy E has a minimum (saddle point), with A = A^ -*, for (t), L, si^) G 7. 
Now, if E is minimal the stability follows by the Liapunov stability theorem ||12|| . If 
E has a saddle point for (ff,L, kip) € 7, and with A = A*- -*, then one argues as fol- 
lows. In any neighbourhood of 7 there exists a solution with E(f], L, kip) < E(^y) and 
A = A^. There are two possibilities: (a) there is a 7^ G {7+, 7-, 71, 72} — {7} such that 
$<(?), L, si >2 ) — > 7^. As any two trajectories have a finite distance, 7 is Liapunov un- 
stable; (b) there exists a series tk — > t + for which g n (§ t (f], L, kip)) — > 0. As <? n (7) = mg, 
by proposition 1, this implies instability. □ 

4.4 Stability and instability for A \JmgI\ 

The theorem of the preceding section, together with the propositions 2-5, completely 
solves the stability problem. However, the answers and criteria are somewhat intricate 
and not very transparent at first sight. Furthermore, in practice, the top will usually be 
launched with an initial rotational energy large as compared to the potential energy, i.e. 
with a value of the conserved quantity A large as compared to -Jmgl\. In this limit the 
stability criteria simplify considerably. We summarize them, for all possible choices of 
the moments of inertia, as they follow from the theorem above. 

Criteria for stability for large A: We distinguish three cases corresponding to three pos- 
sible choices of the moments of inertia: 

(I) I\ > 1^(1 + a): In this case there exists a positive number cq such that for all |A| > Co 
the rotating solutions 7+ and 7_ are Liapunov unstable, and there is exactly one equiv- 
alence class of tumbling motions which is Liapunov stable. 

(II) I\ < 1^(1 —a): In this case there exists a positive number cq such that for all |A| > Cq 
both 7+ and 7_ are stable. There is exactly one equivalence class of tumbling motions 
which, however, is unstable. 

(III) J 3 (l — a) < Ii < 1^(1 + a): In this case there exists a positive number c\ such 
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that for all |A| > c\ the non-inverted, rotating solution 7+ is Liapunov unstable, while 
the completely inverted rotating solution 7_ is Liapunov stable. There are no tumbling 
motions. 

These statements are easily verified: Case (I) follows from parts (i) and (ii) of the theo- 
rem. Indeed, for large |A| the criterion (i) of proposition 3, eq. (39), simplifies to 

7 3 (1 =f a) — I\ > (stability), or < (instability) for y± . 

With < a < 1 we have I\ > I 3 (l + a) > I 3 (l — a). Thus, by (i) 7+ is unstable, and, 
by (ii), 7_ is unstable too. Furthermore, using eq. (28) for 7] 3 in the limit of large |A| the 
expression (40) simplifies to J 3 (l — a 2 ) — I\. With J 3 (l + a) < I\ as assumed we also 
have the inequality J 3 (l — a 2 ) < I\. By (T4) (i) of sec. 3.3 there exists one tumbling 
solution 71 which, by part (iii) of the theorem, is Liapunov stable. 
Case (II) is completely analogous to case (I) but this time 7+ is stable, and so is 7_. 
The tumbling solution 71 (the only one that exists) is unstable. 

In case (III), finally, part (i) of the theorem implies that 7+ is unstable, while part (ii) 
implies that 7_ is stable. Furthermore, the result (T4) (ii) of sec. 3.3 tells us that here 
there is no tumbling solution. This completes the proof of the criteria in the three cases. 
Clearly, case (III) is the genuine tippe top whose strange behaviour we described in the 
introduction and which triggered this analysis. 

5. Numerical results and summary 



In order to illustrate our results we have studied numerical solutions of the equations 
of motion (14) for the three characteristic situations described in sec. 4.4 above [|T0 



It is useful to write the variables L and Si^ in dimensionless form. Having chosen the 
radius r of the sphere to be unity this is achieved by expressing L in units of m^/g, s li2 
as well as any other velocity in units of y/g, time in units of g^ 1 ^ 2 . In what follows all 
variables are given in these rational units, viz. 

L = 1 L u = 1 ■ t = g^_. 

The values of the original physical quantities are recovered by multiplying angular mo- 
menta by m(yf 1 / 2 r 3 / 2 , velocities by (gr) 1 ^ 2 , times by r 1 / 2 ^ 1 / 2 . The relevant constants 
which determine the behaviour of the top are the distance a of the center-of-mass from 
the geometric center, the asymmetry e := (J 3 — of the moments of inertia, and 

c := 1/ii, the inverse of the transversal moment of inertia I2), expressed in units 
of mr 2 . In all examples we have chosen the coefficient to be fi = 0.75, in order to have 
the solutions approach their asymptotics rapidly. 

Figs. 3 and 4 pertain to case (III) of the criteria in sec. 4.4, with a = 0.1, e = 0, 
c = 2.5. Fig. 3 shows the position of the top's symmetry function of time, for 
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an initial condition close to rotation about the (positive) vertical. The rapid oscillations 
which are superimposed to the inversion of the top are easy to understand: they are a 
remnant of the nutational nodding that the top would perform if it were force- free ||. 
Indeed, for a = and fi = the equations of motion (14) simplify to 

dr) 1 „ dL 

whose solution describes uniform rotation of the symmetry axis about the constant 
angular momentum. As a and fi are small in the chosen example, by a theorem of 



Poincare |12] which guarantees a smooth transition of the solution to the force-free 
solution, as a and \x tend to zero, the behaviour of the top still reflects that nutation. 
Fig. 4 shows the time evolution of t)(£), the motion starting at the top of the figure and 
ending in the completely inverted position. 

Figs. 5 and 6 illustrate case (I) of the criteria, i.e. l\ > + a), or e < —a for an 
initial rotation close to the positive vertical (fig. 5), and an initial rotation close to the 
negative vertical (fig. 6). The constants are chosen to be a = 0.1 and e = —0.3. Both 
solutions move quickly towards the tumbling motion 71 which is asymptotically stable. 

The case (II) of the criteria, sec. 4.4, describes an "indifferent" top for which both 
7+ and 7_ are asymptotically stable. This is illustrated by figs. 7 and 8 where we have 
chosen a = 0.1 and e = 0.2. In fig. 7 the top starts at the bottom of the figure and 
moves towards the upright position, the initial value being chosen in the attractive basin 
of 7 + . In fig. 8 the top starts at the top and quickly tends to the inverted position 7_. 

Fig. 9, finally, also belongs to the case (II) which also admits one equivalence class of 
tumbling motions. Unlike 7+ and 7_, however, these are Liapunov unstable. 

In summary, the equations of motion for the tippe top, subject to gravitation and 
to sliding friction, can be formulated in terms of an optimally adapted, minimal set of 
coordinates. The conservation law (2), called " Jelett's integral" in the early literature on 
this topic, and which can be derived by a purely geometric argument ||, follows from 
these equations in an elementary and transparent manner. The total energy is found 
to be a suitable Liapunov function for the stability analysis of the spinning tippe top. 
Its extrema on the hypersurfaces defined by the conservation law A = A^ -* = const, 
are studied. The solutions of constant energy which are the asymptotic states of the 
top, in case of stability, are obtained explicitly from the equations of motion in the 
limit of vanishing sliding friction. Our main result is contained in the theorem of sec. 
4.3 which answers the question of asymptotic Liapunov stability for all choices of the 
constants. The criteria provided by the stability theorem simplify somewhat in case the 
rotational kinetic energy is large. The results are given in sec. 4.4. Finally, numerical 
sample calculations, illustrated by figs. 3-9, confirm the salient features of our analysis. 
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Appendix 



The coefficient g n describing the normal force F n = g n e 3 is calculated as follows. 
One calculates first the acceleration s 3 of the center-of-mass in the vertical direction by 
taking the orbital derivative of eq. (12). The result is 



a 



h = -—g n (r),L,s h2 ) 

(t) x e 3 | [afi x e 3 - ^(afj - e 3 ) x v}} 



~(\Lx[Lxf,\). 



This must be equal, by Newton's law, to g n jm — g. Working this out yields the desired 
formula for g n , 

r . . mgh{l + a{ m L 2 -L,L,)/{gll)} 

g n (rj, L, si, 2 ) = — w £— T t-^ rr-j — - A.l 

h + ma 2 (l - 7?|) + mafi{(r] 3 - a)e 3 - (1 - ar] 3 )r]\ ■ v 
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Figure captions 



Fig. 1 Axially symmetric top of spherical shape (radius r = 1). The center-of-mass S 
is at distance a from the center M. Sliding friction related to rotation about the 
3- or the 3-axis, active in the point of support A, is perpendicular to the plane 
of the drawing. The moment arms of the corresponding torques R and R are as 
indicated. 

Fig. 2 The symmetry axis of the top is described by f), L is the angular momentum, 
both with respect to the system K attached to the center-of-mass whoses axes are 
parallel to the ones of the inertial system K (not shown). L need not be in the 
plane of v) and e 3 . Also not shown are the velocities of S and A w. r. t. K . 

Fig. 3 Type (III) top: time evolution of rj 3 for initial condition i) = (0.4,0, V0.84), L = 
(—1, 0, 5), u(0) = 0. The oscillations show the nutational nodding of the top. 

Fig. 4 Inversion of the tippe top (type (III)) with initial conditions t) = (0.2,0, V0.96), 
L = (0,0,5), u(0) = 0. The top starts near rj = e 3 and moves quickly towards 
V = -e 3 . 

Fig. 5 Motion of a type (I) top towards a tumbling solution, the initial conditions being 
f] = (0.2,0, V / 0T96), L = (0,0,5), u(0) = 0. The top starts near r) = e 3 but 
stabilizes in a tumbling motion. 

Fig. 6 Same case as in fig. 5, except that the top is launched near r) = — e 3 , i.e. with 
initial conditions r) = (0.2, 0, -VOM), L = (0, 0, 5), u(0) = 0. 

Fig. 7 Indifferent top (case (II)) for which both the upright position and the inverted 
position are Liapunov stable. The initial conditions are f) = (0.8,0,0.6), L = 
(0,0,5), u(0) = 0. It moves towards the upright position. 

Fig. 8 Same top as in fig. 7 but launched with initial conditions f] = (V0.84, 0, 0.4), 
L = (0,0,5), u(0) = 0, i.e. in the basin of attraction of the inverted asymptotic 
state. 

Fig. 9 Same top as in figs. 7 and 8, launched with initial conditions 

f] = (VI - 0.495 2 , 0,0.495), L = (0,0,5), u(0) = 0. This state remains in the 
tumbling regime which, however, is Liapunov unstable. 
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